HU-EP-03/54 
September 2003 



m 
o 
o 

(N 

CD 
00 

(N 



The Aoki phase for Nf = 2 Wilson fermions revisited * 
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We report on a numerical reinvestigation of the Aoki phase in full lattice QCD with two flavors of unimproved 
Wilson fermions. For zero temperature the Aoki phase can be confirmed at inverse coupling /3 = 4.0 and /3 — 4.3, 
but not at P — 4.6 and f3 = 5.0. At non-zero temperature the Aoki phase was found to exist also a,t f3 — 4.6. 
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1. The Aoki phase 

In studies of chiral symmetry breaking it is de- 
sirable to start from an (almost) chirally invari- 
ant formulation of lattice fermions. At present 
this is best realized by using fermions satisfy- 
ing the Ginsparg- Wilson relation, e.g. using the 
overlap operator, or by working with the domain 
wall framework [1]. In such approaches Wilson 
fermions are still important as an input. 

For the Wilson-Dirac operator (which breaks 
chiral invariance explicitly) in Nf = 2 lattice 
QCD Aoki [2] has predicted that in a certain pa- 
rameter range there is a phase in which parity 
and flavor symmetry are both spontaneously bro- 
ken. It is separated from an unbroken phase by a 
line of second order phase transition on which the 
pion states are expected to become massless [3]. 
In Fig. 1 the proposed phase diagram is sketched 
in the (/?, K)-plane . 

In a recent paper [4] it has been pointed out 
that the features of the Aoki phase are of rel- 
evance for locality and for restoration of chi- 
ral invariance in quenched and full QCD with 
Ginsparg- Wilson and domain wall fermions. Ac- 
cordingly, the region of the Aoki phase has to be 
avoided in such investigations in order not to spoil 
physical reliability. 

In the light of this it is important to have pre- 
cise information about the region where the Aoki 
phase really occurs. In the past it has been ques- 
tioned [5] whether the Aoki phase exists also at 



larger /3- values, i.e. towards the continuum. How- 
ever, previous investigations of this issue did not 
lead to a unique answer. Therefore we have ad- 
dressed this problem again [6]. 

Here we report on results showing that most 
likely the Aoki phase does not exist for /3 > 4.6 
at zero temperature. Furthermore, we present 
data for finite temperature indicating such bro- 
ken phase still at /3 = 4.6 and connecting it to 
the finite temperature phase transition. 




*Talk presented by A. Sternbeck at the International Sym- 
posium Lattice 2003, Tsukuba, Japan. 



Figure 1. Phase structure proposed by Aoki 
et ai. The shaded region B is the phase where 
flavor and parity are spontaneously broken. Both 
symmetries are conserved in regions A. 



2. Simulation details 

To investigate the Aoki phase we have simu- 
lated lattice QCD with two flavors of unimproved 
Wilson fermions using the Hybrid Monte Carlo al- 
gorithm. An explicitly symmetry-breaking source 
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Figure 2. {■ipij^T^ tp) for f3 = 4.3 as a function of 
K at several values of h (left, the lines are spline 
interpolations to guide the eye) and as a func- 
tion of ft, at K = 0.21 (right, the line represents a 
fit described in the text). The arrow marks the 
extrapolated value of the order parameter. 



0.04 



0.03 



0.02 




0.1 



0.08 



0.06 



0.04 



0.02 



E /3 = 4.6 



0.197 



0.198 



0.199 



6")! = 0.1986 I — < 
8*11 = 0.1984 I — 1 
10* K = 0.1981 'r--< 



0.01 0.02 
h 



0.03 



Figure 3. {ipijsT^i;) for /? = 4.6 at ft = 0.005 
as a function of k, for different lattice sizes (left, 
the lines are drawn to guide the eye), and as a 
function of ft for k changing from 0.1986 to 0.1981 
(right, the line represents a fit described in the 
text). The arrow marks the extrapolated value of 
the order parameter. 



term was added to the fermion matrix Mw, 

M{h)^Mw + hij5T\ (1) 

since without the added term the order parameter 
{'ipij^T^ip) would always be zero on a finite lat- 
tice. Then (Tpij^r^ip) was measured varying the 
lattice size V and the (non-zero) ft-values. The 
order parameter {ipij5T^ip)h=o monitoring spon- 
taneous symmetry breaking is obtained by taking 
the double limit in the following order 

(-0i75r^V)/i=o = lim lim (7/^i75T^'(/') ■ (2) 

For zero temperature the simulations were per- 
formed on lattices ranging from 4'* to 12'* at f3- 
values 4.0, 4.3, 4.6, and 5.0, with k and ft in the 
intervals 0.15 < k < 0.28 and 0.003 < ft < 0.04, 
respectively. For each (3 and each lattice size 
(starting on 6^) the physically relevant values of 
K were determined first, namely the region where 
(Tpij^T^ip) has a peak at finite ft. 

This is shown on the left hand sides of Figs. 2 
and 3. One sees that for /3 = 4.3 and (3 — 4.6 the 
crucial values are near k = 0.21 and k = 0.1981, 



respectively. At these values we have increased 
the lattice sizes until the measurements for dif- 
ferent sizes agreed within errors. In this way we 
achieved that we can treat our largest lattices as 
infinitely large. 

3. The extrapolation h — > 

As described in Ref. [6] a reasonable ansatz to 
extrapolate the data to ft = is 

f{h)^A + Bh^' + ... . (3) 

With C as a free parameter it describes the data 
well. Indeed, the parameter of interest. A, is 
then robust against the introduction of correction 
terms linear and quadratic in ft. 

At (3 = 4.3 the parameter {-ipij^T^ip) has a non- 
vanishing extrapolation, while at (3 = 4.6 the best 
fit gives A = 0. Combined with similar simula- 
tion results at /3 = 4.0 and P = 5.0 (see [6]) we 
confirm the Aoki phase at {f3, k) = (4.0, 0.22) and 
(4.3,0.21). The fit parameters B and C turn out 
to agree for different f3 within errors (the values 
are B « 1 and C « 0.65). At /? = 4.6 the ansatz 
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The area of the phase diagram 
Ref. [6]. Squares indicate that 
{i^i^bT^ tp) h=o is finite, diamonds indicate that 
{'tpi'y5T^4')h=o vanisches, and stars denote points 
where {ipi^^T^ 'ij}) h=ti 7^ is uncertain. The fines 
sketch the hmits of the Aoki phase. The data in- 
dicate that the point (/?, k) = (4.0, 0.215) (circle) 
is very close to the lower critical line [6]. 



also works, but results in a vanishing order pa- 
rameter. Also at /3 = 5.0 no finite order param- 
eter was found [6]. This suggest that the Aoki 
phase ends somewhere near (/3, n) = (4.6, 0.1981) 
as shown in Fig. 4. The shift in k with increasing 
lattice size (see Fig. 3) appears to be related to 
the end of the symmetry-breaking phase. 

4. Outlook for finite temperature 

We have also performed simulations for the 
finite-temperature case. In contrast to zero tem- 
perature we find the Aoki phase aX (3 — 4.6 
as shown in Fig. 5. One can see that the 
Polyakov loop (open symbols) steeply rises near 
K — 0.19705 where {■ipij^T^^') (filled symbols) ex- 
trapolates to a finite value. This means that the 
finite-temperature phase transition occurs close 
to the Aoki phase. We were not able to detect a 
separation of this transition from the Aoki phase. 
This is in contrast to Ref. [7] where a scenario 
with separated transition lines is proposed. 
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Figure 5. The parameter {ipij^r^ip) ^^'^ the 
Polyakov loop as function of k on a 8'^ x 4 and a 
10^ X 4 lattice at ft. = 0.005 and P = 4.6 (left, the 
lines are drawn to guide the eye). At k = 0.19705 
we obtained {'tpi'y^T^il;) ^ from an extrapolation 
explained in the text (right). 
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